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Abstract. We look for space-time correlations between 
2233 gamma-bursts in the Batse Catalogs and 894 
upward-going muons produced by neutrino interactions 
in the rock below or inside MACRO. Considering a search 
cone of 10° around GRB directions and a time window 
of ±200 s we find events to be compared to 0.035 ex- 
pected background events due to atmospheric neutrinos. 
The corresponding upper limit (90% c.l.) is 0.87 x 10 -9 
cm~ 2 upward-going muons per average burst. 
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1. Neutrino astronomy with MACRO 

Besides high energy gamma ray production resulting from 
7T° decay, astrophysical beam dump models predict neu- 
, trino emission from decay. Mesons are produced by 
accelerated protons interacting with matter or photons in 
an accretion disk (Gaisser 1995] ). The discovery of TeV 
gamma-ray emissions has enhanced the potential possi- 
bilities of this mechanism and the possible existence of 
such sources, but energies are not high enough to es- 
clude syncroton radiation or bremsstrahlung and inverse 
Compton production mechanisms. Neglecting photon ab- 
sorbtion, it is expected that neutrino fluxes are almost 
equal to gamma ray ones and that the spectrum has the 
typical form due to the Fermi acceleration mechanism: 
4j| oc E~ 2 - ^ 2 - 2 . Neutrinos produced in atmospheric cas- 
cades are background to the search for astrophysical neu- 
trinos for energies < 10 TeV. In fact, atmospheric neutri- 
nos have a softer spectrum than astrophysical neutrinos, 
since at energies > 100 GeV the decay length of mesons 
in the atmosphere becomes longer than atmospheric depth 
and the spectrum steepens (differential spectral index 7 ~ 
3.7). GRBs are possible sources of high energy vs: in the 
fireball scenario the beam dump mechanism can lead to v 
emission (Halzen 1996, Waxman 1997). The trancience of 
the GRB emission, even though constrains the maximum 



neutrino energy to < 10 19 eV, improves GRB association 
with measured neutrino events in underground detectors 
using arrival directions and times, even though expected 
fluxes are much lower than atmospheric v background. 

The MACRO detector, located in the Hall B of the 
Gran Sasso underground laboratories, with a surface of 
76.6 x 12 m 2 and a height of 9 m, can indirectly detect 
neutrinos using a system of ~ 600 ton of liquid scintilla- 
tor to measure the time of flight of particles (resolution 
~ 500 ps) and ~ 20000 m 2 of streamer tubes for track- 
ing (angular resolution better than 1° and pointing ac- 
curacy checked using moon shadow detection (Ambrosio 



1998a )). The time of flight technique allows the discrim- 



ination between downward-going atmospheric muons and 
upward-going events produced in the rock below (aver- 
age atmospheric neutrino energy (E u ) ~ 100 GeV) and 
inside ((-Ev) ~ 4 GeV) the detector by neutrinos which 
have crossed the Earth. Between ~ 31 x 10 6 atmospheric 
muons, a sample of 909 upward-going muons is selected 
with an automated analysis. The data taking has begun 
since March 1989 with the incomplete detector (Ahlen 
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1995) and since April 1994 with the full detector (Ambro- 
sio 19981). In our convention 1/(3 = ATc/L, calculated 
from the measured time of flight At and the track length 
between the scintillator layers, is ~ 1 for downward-going 
muons and ~ — 1 for upward-going muons. Events with 
-1.25 < l/P < -0.75 are selected. 

We look for a statistically significant excess of v events 
in the direction of known 7 and AT-ray sources (a list of 
40 selected sources, 129 sources of the 2 nd Egret Cata- 
logue, 2233 Batse GRB, 220 SN remnants, 7 sources with 
7 emission above 1 TeV) with respect to fluctuations of 
the atmospheric v background. The expected background 
from atmospheric vs is calculated in declination bands of 
±5° around the declination of the source mixing for 100 
times local coordinates and times of upward-going jtts. We 
calculate flux limits in half-cones of 3° taking into account 
reduction factors due to the signal fraction lost outside 
the cone (which depends on v spectrum, kinematics of 
CC interaction, /i propagation in the rock, MACRO an- 
gular resolution) . We do not find any signal evidence from 
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known sources or of clustering of events (we measure 89 
clusters of > 3 events and expect 81.2 of them in a cone 
of 3°). Muon flux limits for some sources are: 2.5 for Crab 
Nebula, 5.6 for MRK421, 3.71 for Her X-l, 0.45 for Vela 
Pulsar, 0.65 for SN1006 in units of xl0~ 14 cm~ 2 s -1 . For 
most of the considered sources MACRO gives the best flux 
limits compared to other underground experiments. 

2. Space-time correlations between GRBs and 
upward-going muons 

We look for correlations with 2233 GRBs in the Batse 



Catalogs 3B and 4B (Meegan |1997| ) collected since 21 
Apr. 1991 to 5 Oct. 1998 and 894 of the 909 upward- 
going muons detected by MACRO during this period (see 
Fig. Q). Considering Batse angular accuracy, we estimate 
that a half-cone of 20° (10°) contains 99.8% (96.8%) of 
v sources (if GRB sources are v sources, too). The same 
percentages of upward-going muons are contained in these 
half-cones from the GRB sources. As a matter of fact, 
we calculate via Monte Carlo the fraction of signal lost, 
which depends on the v spectral index, multiple scatter- 
ing of muons during propagation in the rock and MACRO 
angular resolution. Using various cone apertures, we esti- 
mate that this fraction is negligible for > 10°. 

The area for upgoing muon detection in the direction of 
the GRBs averaged over all the bursts is 119 m 2 . Its value 
is small because MACRO is sensitive to neutrinos only in 
the lower hemisphere and because it was incomplete in the 
period 1991-1994. We find no statistically significant cor- 
relation between neutrino event and GRB directions and 
detection times. As shown in Fig. ^, we find no events in 
a window of ±200 inside 10° from GRB directions and 1 
event inside 20°, which was measured after 39.4 s from the 
Batse GRB of 22 Sep. 1995 (4B 950922). For this burst 
the radius of the positional error box in the Batse catalog 
is 3.86°, much smaller than the angular distance of 17.6° 
at which we find the neutrino event. The expected number 
of atmospheric v events is computed with the delayed co- 
incidence technique. We expect 0.035 (0.075) events in 10° 
(20°). The corresponding upper limits (90% c.l.) for the 
upward-going muon flux are 0.87 x 10~ 9 cm~ 2 (1.44 x 10~ 9 
cm~ 2 ) upward-going muons per average burst. These lim- 
its exclude an extreme cosmic string-type model reported 



in (Halzen 1996), which results in 10 1 /i cm 2 , while ac- 



cording to a fireball scenario model in (Waxman 1997), a 
burst at a distance of 100 Mpc producing 0.4 x 10 51 erg in 
neutrinos of around 10 14 eV would produce ~ 6 x 10~ n 
cm -2 upward-going muons. 
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Fig. 1. MACRO upward-going muon (dashed line) and 
Batse GRB (solid line) distributions versus year. 
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Fig. 2. Difference in detection times vs cosine of angular 
separation between Batse GRBs and upward-going /is. (a) 
and (b) have different time scales. In (a) the ±200 s - 
10°, 20° windows are indicated. 



Ambrosio, M., et al., MACRO Coll., 1008, Phys. Lett. B434 
451 

Gaisser, T.K., et al, 1995, Phys. Rep. 258, 173. 
Halzen, F., Jaczko, G., et al, 1906, Phys. Rev D54, 2779 
Meegan ; C.A., 1997, Current Batse Gamma-Ray Burst C ata- 
log, http://www.batse.msfc.nasa.gov/data.grb/cataloj 



Waxman, E., Bahcall, J., 1907, Phys. Rev. Lett. 78, 2202; hep- 



shaitinwing of rns mir rays by thp M non using a dppp under- 



ph/0807282 



ground detector, hep-ex/0807006, accepted by Phys. Rev 



D 

Ahlen, S., et al., MACRO Coll., 1995, Phys. Lett. B357 481 



(*) The MACRO Collaboration 



Montaruli T. & Ronga F.: Space-time correlation between vs and GRBs 



3 



M. Ambrosio 12 , R. Antolini 7 , C. Aramo 7 '™, 
G. Auricmma 14 ' a , A. Baldini 13 , G. C. Barbarino 12 , 

B. C. Barish 4 , G. Battistoni 6 ' 6 , R. Bcllotti 1 , 

C. Bemporad 13 , E. Bcrnardini 2 ' 7 , P. Bcrnardini 10 , 
H. Bilokon 6 , V. Bisi 16 , C. Bloise 6 , C. Bower 8 , 
S. Bussino 14 , F. Cafagna 1 , M. Calicchio 1 , D. Campana 
12 , M. Carboni 6 , M. Castellano 1 , S. Cecchini 2 < c , F. Cei 
n < 13 , V. Chiarella 6 , B. C. Choudhary 4 , S. Coutu n <°, 
L. De Benedicts 1 , G. De Cataldo 1 , H. Dekhissi 2 ' 17 , 

C. De Marzo 1 , I. De Mitri 9 , J. Derkaoui 2 < 17 , 
M. De Vincenzi 14 ' 6 , A. Di Credico 7 , O. Erriqucz 1 , 
C. Favuzzi 1 , C. Forti 6 , P. Fusco 1 , G. Giacomelli 2 , 

G. Giannini 13 ^, N. Giglictto 1 , M. Giorgini 2 , M. Grassi 
13 , L. Gray 7 , A. Grillo 7 , F. Guarino 12 , P. Guarnaccia 1 , 

C. Gustavino 7 , A. Habig 3 , K. Hanson n , R. Heinz 8 , 
Y. Huang 4 , E. Iarocci 6,9 , E. Katsavounidis 4 , E. Kearns 
3 , H. Kim 4 , S. Kyriazopoulou 4 , E. Lamanna 14 , C. Lane 

5 , D. S. Levin 11 , P. Lipari 14 , N. P. Longley 4 < ( , 
M. J. Longo n , F. Maaroufi 2 ' 17 , G. Mancarella 10 , 
G. Mandrioli 2 , S. Manzoor 2 < m , A. Margiotta Neri 2 , 
A. Marini 6 , D. Martello 10 , A. Marzari-Chiesa 16 , 
M. N. Mazziotta 1 , C. Mazzotta 10 , D. G. Michael 4 , 
S. Mikheyev 4 > 7 <'\ L. Miller 8 , P. Monacelli 9 , 
T. Montaruli 1 , M. Monteno 16 , S. Mufson 8 , J. Musser 8 , 

D. Nicol6 13 ' d , C. Orth 3 , G. Osteria 12 , M. Ouchrif 2 < 17 , 
O. Palamara 10 , V. Patera 6 ' 9 , L. Patrizii 2 , R. Pazzi 13 , 

C. W. Peck 4 , S. Petrera 9 , P. Pistilli 14 < e , V. Popa 2 <\ 
A. Raino 1 , A. Rastelli 2 ' 7 , J. Reynoldson 7 , F. Ronga 6 , 
C. Satriano 14 ' a , L. Satta 6 ' 9 , E. Scapparone 7 , 
K. Scholberg 3 , A. Sciubba 6,s , P. Serra-Lugaresi 2 , 
M. Severi 14 , M. Sioli 2 , M. Sitta 16 , P. Spinelli 1 , 
M. Spinetti 6 , M. Spurio 2 , R. Steinberg 5 , J. L. Stone 3 , 
L. R. Sulak 3 , A. Surdo 10 , G. Tarle 11 , V. Togo 2 , 
D. Ugolotti 2 , M. Vakili 15 , C. W. Walter 3 , and 
R. Webb 15 . 



12. Dipartimento di Fisica dell'Universita di Napoli and INFN, 
80125 Napoli, Italy 

13. Dipartimento di Fisica dell'Universita di Pisa and INFN, 
56010 Pisa, Italy 

14. Dipartimento di Fisica dell'Universita di Roma "La 
Sapienza" and INFN, 00185 Roma, Italy 

15. Physics Department, Texas A&M University, College Sta- 
tion, TX 77843, USA 

16. Dipartimento di Fisica Sperimentale dell'Universita di 
Torino and INFN, 10125 Torino, Italy 

17. L.P.T.P., Faculty of Sciences, University Mohamed I, B.P. 
524 Oujda, Morocco 

a Also Universita della Basilicata, 85100 Potenza, Italy 

b Also INFN Milano, 20133 Milano, Italy 

c Also Istituto TESRE/CNR, 40129 Bologna, Italy 

d Also Scuola Normale Superiore di Pisa, 56010 Pisa, Italy 

e Also Dipartimento di Fisica, Universita di Roma Tre, Roma, 

Italy 

/ Also Universita di Trieste and INFN, 34100 Trieste, Italy 
g Also Dipartimento di Energetica, Universita di Roma, 00185 
Roma, Italy 

h Also Institute for Nuclear Research, Russian Academy of Sci- 
ence, 117312 Moscow, Russia 

i Also Institute for Space Sciences, 76900 Bucharest, Romania 

I Swarthmore College, Swarthmore, PA 19081, USA 

m RPD, PINSTECH, P.O. Nilore, Islamabad, Pakistan 

n Also INFN Catania, 95129 Catania, Italy 

o Also Department of Physics, Pennsylvania State University, 

University Park, PA 16801, USA 



1. Dipartimento di Fisica dell'Universita di Bari and INFN, 
70126 Bari, Italy 

2. Dipartimento di Fisica dell'Universita di Bologna and INFN, 
40126 Bologna, Italy 

3. Physics Department, Boston University, Boston, MA 02215, 
USA 

4. California Institute of Technology, Pasadena, CA 91125, 
USA 

5. Department of Physics, Drexel University, Philadelphia, PA 
19104, USA 

6. Laboratori Nazionali di Frascati dell'INFN, 00044 Frascati 
(Roma), Italy 

7. Laboratori Nazionali del Gran Sasso dell'INFN, 67010 As- 
sergi (L'Aquila), Italy 

8. Depts. of Physics and of Astronomy, Indiana University, 
Bloomington, IN 47405, USA 

9. Dipartimento di Fisica dell'Universita delPAquila and INFN, 
67100 L'Aquila, Italy 

10. Dipartimento di Fisica dell'Universita di Lecce and INFN, 
73100 Lecce, Italy 

11. Department of Physics, University of Michigan, Ann Ar- 
bor, MI 48109, USA 



